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Abstract 

We compute the correlation functions mixing the powers of two non-commuting ran- 
dom matrices within the same trace. The angular part of the integration was partially 
known in the literature |16II17| : we pursue the calculation and carry out the eigenvalue 
integration reducing the problem to the construction of the associated biorthogonal 
polynomials. The generating function of these correlations becomes then a determi- 
nant involving the recursion coefficients of the biorthogonal polynomials. 



1 Introduction and main result 

Random matrix models were first introduced in the context of nuclear physics in order to describe the energy 
levels statistics for very large nuclei. Wigner proposed the hypothesis that these were distributed as the 
eigenvalues of a matrix with random entries. Later random matrix models were used in many areas of physics 
and mathematics [HI H2| • 

An important application of random matrices is to 2d gravity, that is, statistical physics on a random surface. 
In fact, the perturbative expansion of a matrix integral can be accomplished by drawing Feynman graphs on 
fixed-genus surfaces. Matrix integration can therefore encode the summation over the set of discretized surfaces 
(possibly carrying some type of matter). 

When the parameters of the model are fine-tuned near a critical point the average graph's size diverges 
and macroscopic graphs dominate the sum, so that suitable critical limits can represent statistical models over 
smooth surfaces. The smooth surfaces that one describes with the aid of matrix models have properties of 
scale invariance: this means that the critical points of matrix integrals are related to representations of the 
conformal group. We recall that its finite dimensional representations are classified by two integers (p,q): it 
is known that one-matrix-models provide instances of (p, 2)-irreps only, where as two-matrix-models allow to 
obtain (p, (^-representations. 
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Possibly the first such application of the two-matrix-model was to describe the Ising model on a random 
surface; in this case the Ising ferromagnetic transition corresponds to the conformal minimal model (3,4). To 
see this, one should associate a color (or spin) + or — to each matrix so that the vertices of one matrix are 
labelled with a plus sign and the vertices of the other matrix with a minus sign. Then the Feynman graphs 
generated by a two matrix model represent discrete surfaces carrying spins (+ and -), i.e. an Ising model on a 
random surface. 

The correlation functions of random matrices generate discrete surfaces with boundaries, and thus are in 
relationship with boundary conformal theory. A formula for correlation functions representing surfaces with 
mono-colored boundaries has been known since JOJ- It is the aim of this paper to give a formula for a mixed 
correlation function, i.e. the generating function for discrete surfaces with a bi-colored boundary. 

The 2-matrix model has attracted a lot of attention recently, and important progress have been made in 
the study of the associated bi-orthogonal polynomials [31 03 EH- Here, we will express the mixed correlation 
function in terms of the bi-orthogonal polynomials. 



1.1 Definition and notation 

We consider two N x N Hermitian matrices Mi , M 2 , with a probability measure 



dn(M 1 ,M 2 ) := Z N ~ 1 dM 1 dM 2 exp 



TV(Vi(Mi) + V 2 (M 2 ) - M X M 2 ) 



dMidM 2 exp 



Tr(y x (Mi) + V 2 (M 2 ) - M X M 2 ) 



;i-i) 



where dM\dM 2 is the product of Lebesgue measures of all the independent real and imaginary parts of the 
components of the two matrices divided by the square of the volume of the unitary group U(N) (for later 
convenience). The functions V\ and V 2 are called the potentials and must be chosen so as to make the integral 
convergent. The normalization factor Zfq is called the "partition function", where the name "function" refers 
to its dependence on the two potentials. 

One can rewrite the measure in term of eigenvalues and angular integrals |12| : 



dfi(M u M 2 ) := ZaT 1 A(X) 2 A(r) 2 exp 



N 



5>ite) + ^(y<)) 



i=l 



N 



M^XUV^YV 



dUdVUdxidyi , 



i=l 



X := diag(xi, . . . ,x N ), Y := diag(y x , . . . , y N ), U,VeU(N) 

r r N r N 

Z N-=J J Yldx^&iXfAiYfdUdVexp - ^(Fi(^) + V 2 (Vi)) 



JJ^XUV^YV 



(1-2) 



where dUdV is the product of the normalized Haar measure over U(N) x U(N). In the original two Hermitian 
matrix model, the integration path for the Xj's and y^s is the real axis (location of the eigenvalues of a Hermitian 
matrix), however, the model can be generalized to include complex paths or their homology classes in case the 
potentials are holomorphic or meromorphic [UEI]. 



1.2 Correlation functions 

In applications of the two-matrix-model to statistical physics on a random surface, one is interested in comput- 
ing correlation functions involving traces of products of powers of M\ and M 2 . Each such correlation function 
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can be expanded in Feynman graphs, which represent discrete surfaces with boundaries: the number of bound- 
aries is the number of traces, the length of each boundary is the total power of Mi plus the total power of M2 
within the trace (Hj- 

For instance: (TrMf) is the generating function for discrete surfaces with one boundary (a circle) of length 
r made of + spins only. (Tr M|) is the generating function for discrete surfaces with one boundary of length s 
made of — spins only. (TrM^MI) is the generating function for discrete surfaces with one bi-colored boundary 
of length r + s made of r + spins, followed by s — spins. 

More generally, (TtMp M^ 1 Mp 2 M| 2 . . . M^Mi^ 1 ) is the generating function for discrete surfaces with one 
2n-colored boundary of length Y^i v\ + ^ Sj made of r\ + spins, followed by s\ — spins, followed by + 
spins, . . ., followed by s n — spins. 

One may also be interested in "multi-loop" correlators (i.e. more than one boundary), for instance: 
(TrM^TrMp) is the generating function for discrete surfaces with two spin + boundaries, one of length r%, 
the other of length r<i- More generally, one may consider correlation functions involving an arbitrary number 
of traces, each containing arbitrary words of M\ and Mi. 

The correlation functions, with an arbitrary number of traces, with each trace containing powers of only one 
matrix have been known since the work of They can be expressed in terms of bi-orthogonal polynomials. 

The aim of the present article is to express the mixed correlation function 

(Tr Ml Ml) (1-3) 

in terms of bi-orthogonal polynomials too, and confirm that the key property of these models is that all relevant 
spectral statistics can be reduced to the computation of the corresponding biorthogonal polynomials 

1.3 Bi-orthogonal polynomials 

Two sequences of monic polynomials 

7r n (x) = x n + ■ ■ • , a n (y) =y n + -.., n = 0, 1, . . . (1-4) 
are called biorthogonal if they are "orthogonal" with respect to a coupled measure on the product space: 

fdxdyn n (x)a m {y)e- v ^- v *W +xy = h n S mn , h n ^ Vn G N (1-5) 



where V\{x) and Viiy) are the functions (called potentials) appearing in the two-matrix model measure (|l-lj) . 
It is convenient to introduce the associated quasipolynomial differentials defined by the formulas 

^ n (x) := —^^{x^^dx (1-6) 

V h n -\ 

My) ■= -^=o-n-i{y)e- V2{y) dy ■ (1-7) 

V "n-l 

In terms of these two sequences of differentials the multiplications by x and y respectively are represented by 
semiinfinite square matrices Q = [Qij]i,jen* and P = [Pij]ijeN* according to the formulae 

X1p n ( x ) = ^2 Qn,m1pm(x) ! #«(?/) = P ™,n4>m{y) 
m m 

Qn,m = = P m n , if n > m + 1. (1-8) 
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The matrices P and Q have a rich structure and satisfy the "string equation" [P, Q] = id. However we do not 
need any of their properties except for eq. (|l-8[) to derive our present results, and therefore we refer for further 
details to PI El El E] where these models are studied especially in the case of polynomial potentials. We also 
point out that the model can easily be generalized to accommodate contours of integration other than the real 
axes [21 [3] leaving intact all the properties which are relevant to the following computations. 

1.4 The main result 

Our goal is to prove a formula for the generating function of the correlators 



(Tr(M 1 r M 2 s ))y 1 y 2 := -J- f dM 1 dM 2 Tr(M 1 r M 2 s ) exp (-Tr (Fi(Mi) + V 2 (M 2 ) - MiM 2 )) . 
AW J 

By generating function we mean the formal double Laurent series 

Tr ( \ /T 1 M )) ■= E^V^WM/M/))^- 
\ \ x ~ Miy- M 2 ) I Vl y 2 



(1-9) 



(1-10) 



The main obstacle to this sort of computations so far was posed by the "angular integration" over the unitary 
group U(N). 

One can trace in the literature various attempts at this computation using the loop equations [3 |S1 E] • A 
closed formula was found in the large N limit in |Sj, but an exact formula for finite N was never derived. 

Our strategy is that of reducing the computation of (|1-9|) or even better ()l-l(jj) to the computation of the 
corresponding biorthogonal polynomials associated with the measure (11-1(1 . 

We can now write down the main result of the paper which is both simple and beautiful: 



Tr 



1 



1 



x - Mi y - M 2 



1 - det 



Vi,V 2 



1 N 



7T 



Nx-Qy-PN 



(1-11) 



where P and Q are the matrices in eq. ()1-8|) and vr denotes the projector 



C N onto the span of the 



N 



first N canonical basis vectors, i.e., the N x oo matrix with nonzero entries vr^j = 1, i = 1, . . . N. Formula 

(|1-11|) should be properly understood in the sense of an identity of formal Laurent series in the variables x and 
y, although it would be possible to give an analytic meaning to both sides. For example from Ijl-llJI one can 
easily obtain the following identities 

1 

(Tr{M\ M 2 )) Vu y 2 = Tr(Q^Pn) - - £ (Tr(Q :, 'n)Tr(Q r ' -1 ^n) - Tr(QW^n)) (1-12) 

3=0 



Tr(Q r P 2 n)-E (Tr(Q kl PU)Tr(Q k2 U) -Tr(Q kl PUQ k2 U)^j (1-13) 



fci+fc2=r— 1 



<Tr(M/M 2 2 )> yiy2 

+ V ( -Tr(nQ fci nQ fc2 nQ fc3 )--Tr(nQ fci nQ fc2 )Tr(ng fc3 ) + -Tr(nQ fcl )Tr(nQ fc2 )Tr(nQ fc3 ) ] (1-14) 
^-^ V 3 2 6 / 

where now II = 7r*7r is the semiinfinite square matrix with nonzero entries II4 j = 1, i = 1, . . . N. 
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In fact in the course of the proof of eq. we will also prove the following strong result for the 

correlations of {Tr(Mi r M 2 s ))y 1 y 2 which naturally extends to arbitrary (analytic) functions f(x),g(y): 



{Tr(M 1 r M 2 r )) Vl}V2 



(1-15) 



/ 1 



det 



min(r,s) 

£ 

n=0 



(-1)" 

(ra + 1)! 



X n +1 



„n-l 



\ 1 

\ x\ 



p n-l 
l n+l 

c n+l 



det 



2/1 



yr 1 



1 \ 

y n +i 



/-I 



det 



1 \ 

X n +l 



det 



f 1 

yi 



b n+l 



Vl 



1 \ 

y n +i 

2/n+l ' 



det 



Ki2(x k ,ye) 



(1-16) 

fc/<n+l 



The formula can be extended to arbitrary analytic functions f(x) and g{y) to give 
(Tr(f(M 1 )g(M 2 ))) Vl y 2 = 



(1-17) 



(-1)" 

^ (n + 1)! 

n=0 V 7 



det 


r 1 ... 1\ 

x\ ■■■ x n+1 

rr n ~ X . . . T n_1 
X l X n+1 

^/(^l)' • 7(^n,+l)> 


det 


f 1 ■•• IX 

2/1 ••• 2/n+i 

2/i 2/n+l 

^(yi)- • -g{yn+i)' 


det 


r 1 ^ 

xx ■■■ x n+ x 


det 


/l ••• 1 X 

2/1 • • • 2/n+l 



det 



^n+l 



Ki 2 (x k ,ye) 



(1-18) 

fc,£<n+l 



2/n+l • 



Here we have used the kernel constructed from the quasipolynomials jlUj . 

^ vr ( x ) a ( V ) e -Vi(x)-V 2 (y) 
K 12 (x,y) := = £ *"W g "W e dxdy 



n=l 



n=0 



(1-19) 



2 Proofs of formulae (ITTB and (ITTO 

We want to compute the expectation values 



(Tr(Mx r M 2 s )) v . v :-. y 



N 



d/i(M 1 ,M 2 )Tr(M 1 r M 2 s ). 



(2-1) 



We denote by {xi\i=x..N and {yi}i=x..N the spectra of the two matrices Mi, M 2 and by U € U(N) the relative 
angles. Then we can reduce the integral (|2-lj) to an integral over the spectra of Mx, M 2 and the unitary group 
of the relative angles. Indeed we have 



{Tr(MlMl)) Vl y 2 :-- 



1 f N 

^— J Y]_dxidyi 



1 



Zjsi 

-Vi( Xi )-V2(v 



d/i(Mi,M 2 )Tr(M 1 r M|; 



(2-2) 



N 

^A 2 (X)A 2 (Y) y XiV / dC/l^fe^^^) , (2-3) 

f-^X JU(N) 
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where X = diag(xi, . . . ,xn), Y = diag(yi, . . . ,un) and A(X), A(Y) denote the Vandermonde determinants. 
The computation requires the knowledge of the two-points correlators for the unitary integral in eq. Q2-3J) 
which we analyze in the next subsection. 



2.1 Two— points correlator for the unitary integral 



The computation of this quantity has been considered in the finite N regime in the two papers |16l I17j . In |17| 
was described a complete algorithm that allows the construction of a formula for the most general correlator 



JU(N) 



t Tr(XWYU) 



(2-4) 



Such algorithm involves the introduction of the parametrization of the unitary group given by the Gel'fand- 
Tsetlin coordinates associated to one of the two matrices Mf, however the computation is highly involved and 
in |17j the problem was not carried through to a completely manageable formula. 

On the other hand in ^] a very simple closed formula was proposed for the two-point correlator (|L^i| 2 )c/(Ar) 
in terms of a generating function, that is 



(2-5) 



/ 1 



det 



A(X)A(F) 



E <p) 



N-l 



n=0 ii<i2<...<in+i ' 



1 \ 



x n ~ l 



/ 1 

y P (h) 



det 



.n-l 



det 



1 \ 



y p(h) 
\ b p(h) 



x 



n 



X in + 1 



X *n + 1 / 



det 



/ 1 

Vp(h) 

\ y p(ii) 



1 \ 

yp{in+l) 



n— 1 

yp{i n+1 ) 

b p(in + l) / 



1 \ 

yp(i n+1 ) 

yp(i n+ i) / 



(2-6) 



with the understanding that the Haar measure of the unitary group U(N) has been normalized to unity. This 
formula will be the starting point of our analysis: however the author of ^H] did not actually prove the formula 
but just made an educated (and -as it turns out- correct) guess. 

Therefore, before proceeding to proving our main result (jl-llj) we want to fill in the gaps between the full 
but unpractical algorithm given in and the practical but unproven formula in |16j . We do not need the full 
generality of JZj: our departure point is formula (1.4) ibidem, restricted to the particular case of the two-point 
correlator. For the ease of the reader we rewrite the aforementioned formula in the notation of our present 
paper 



N-l 



N-l 



(Ui h ul l .--u ljn ulj v{N) 



Oq h 

Jn ">n 



A(X)A(y 2 ,...,y N ) 



'N-l 

n 

\k=l 



1=1 



Y[(x e - x h ) ■ ■ ■ Yl(x e 

&31 tyjn 



N 



N-l 



x exp 



E x ^ ~ E ^ kyi 



det 



ik=l 



k=l 



i,j=l..N-l 



(2-7) 
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Here the Haar measure of the unitary group has been normalized to unity, A(y 2 , • • • ,Un) denotes the Van- 
dermonde determinant of the N — 1 numbers 2/2, •••yjv, while A(X) is the short form for the Vandermonde 
determinant of the whole spectrum of the matrix X. 

Eq. 1)2-7(1 was proven rigorously in for x\ < x 2 < ■ ■ ■ < xtv; however it is straightforward to realize 
that the unitary integral in eq. (|2-3|) defines an analytic (in fact entire) function of the variables X and Y. 
Therefore the result extends to statistical ensembles of pairs Mi, M2 of normal matrices 3 with their spectrum 
on arbitrary paths on the complex plane. Moreover the restriction on the order of the spectrum can be lifted 
because the result is analytic in the variables x\, . . . , x^ and can be analytically continued to C . 

We now set n = 1 in eq. (|2-7|) and hence k\ = j\ = i. Then we have : 



m 1 



3 L x m 



N-l 



u\ W{N) = — — -pr E^II/ d^ k - Xi )e^M^~^= (2-f 

A(X)A(y 2 , . . . , y N ) [[{xe - a*) pes^ k=i Jx * 



3 L BiVl 



A(X)A(y 2 , ...,y N ) i\^(x e - x t ) ^ e ^ 



Xk+l Xi 



y P (kyn-yi {y P {k)H-vif 



Xfc Xi 



y P (k}H-yi {y P {kyri-yif 



A(X)A(Y) n^(x, - ~x~) H^fe - yi) ^ <P) 



(x k+ i ~ Xi)(y p{m - yi) - lVwk,w + i-Wi) _ (( Xk - Xl )(y p(m - yi) - lje°*<VpW-vi) 
We observe that the following identity holds 



(x i -x i )(y m -yi)-Ae xm 



det 



(2-9) 
=(2-10) 
(2-11) 
=(2-12) 

= (2-13) 
(2-14) 



e,m=l...N 



which is realized by performing elementary row operations on the matrix inside the determinant in Q2-14j) so 
that the N x N determinant reduces to a (N — l) 2 determinant by use of Laplace's formula. 

The more general case of the expectation value of the element is obtained by permutation of the 

spectrum of the matrix Y so as to give the formula 



(\ U ji\ 2 )u(N) 



1 



det 



{xt>-Xi){y m -y j )-l)e Xty 



\m=l...N 



A(X)A(Y) 



Y[(x£ - x^ Y[(ye - Vj) 



(2-15) 



This form of Shatashvili's formula (|2-7|l for the case n = 1 is remarkably simple but not suitable for our later 
purposes. Moreover it is not yet clearly equivalent to Morozov's formula (|2-6|) . which is what we need for our 
computation. It will be proved in appendix EI that the two formulae are indeed equivalent. 

3 We recall that a normal matrix is a matrix that commutes with its Hermitian-transposed. Any such matrix can be diagonalized 
using a unitary transformation (and vice- versa). 
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2.2 The correlators (TrMi r M 2 ' 5 )y li v 2 and their generating function 

Starting from eq. (|2-3|) and using formula (|2-6|) with etj = xf and bj = yf we obtain 

N 



[ AUTi{X r U^Y s U)^ xu ^ = V Xi T yf [ dU\Uji 

JU(N) i . j 7t/(JV) 



2 Tr(Xl/tY[/) 



(2-16) 



/ 1 



X j -i 



det 



I \ / 1 
det ' 



n-1 



A(X) 



1 _ _ N ~ 1 \ r 



n-l 



V 4 • • • / 



1 \ 

Hpiin+l) 



..n-l 



P&S N 



n=0 ii<i2<---<j„+i 

det 



1 \ 



^p(ii) 



y 



n-l 

P(*n + l) 



det 



fej ^p(wi) / ,(2-17) 



/I ... 1 \ 

y p{ n) ••• y P ( in+l) 



with the understandings that we use the normalized Haar measure over the unitary group U(N) and that for 
n = the ratio of determinants should be x r , y s ,. •,. The first observation is that the sum over n does not 
actually need to be extended up to the size N of the random matrices because the determinants will vanish 
for n > min(r, s). The next remark is that the ratios of determinants actually define certain totally symmetric 
polynomials of their arguments of degree r — n and s — n respectively: in fact they are Schur polynomials 
corresponding to hook Young diagrams 

/ 1 ••• 1 \ 



det 



SriXix > • • • j x i n +i) : " 



T ™-1 . . . 

X «l X in + 1 



( i 



e n 



det 



1 \ 11-* 

ai<a2<...<a r _ n fc= 1 jiH hjn+l=r-n 



(2-18) 



and a similar expression for the y part. It is interesting to notice that in eq. (|2-17|) the characters of the 
representations of the group GL(N) appear; the same equation could possibly be derived from the character 
expansion of the integrand. 

The formal generating function of these Schur polynomials is: 



oo n+1 ^ 



r=0 



fe=l Z X * fc 



(2-19) 



Eq. Q2-3JI now becomes 



d/z(Mi,M 2 )Tr(MfM|) 



v 



[] d^d^) A(X)A(y)^e( /0 )e^ x 

peS N 



i=l 
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min(r,s) 

x E ( _1 ) n E S r{x h ,...,x in+1 )S s {y p{il) ,...,y p{in+l) ), (2-20) 

n=0 u<i2<...<«n+i 



where d/i(x) := e Vl ^dx and du(y) '■= e ^ 2 ^dy. 
Using eq. (|2-19|) we have: 



J MMi, 



M 2 )Tr 



1 1 



x - Mi y - M 2 , 

AT-1 n+1 



J] d^x^dKm) A(x) A(y £ (-!)« 53 [J (x _ x . )(V _ V( ~) ■ ( 2 - 21 ) 

i=l peSW n=0 ii<i 2 <...<i„+i &=1 * fc i/ / 9 ( 4 fc) ; 

By a relabelling of the y's the sum over p becomes an overcounting factor iV!: 

1 1 



d/x(M!,M 2 )TV 



x - Mi y - M 2/ 

TV AT-1 n+1 



II WMw) A(A:)AGK)eS— £ II (g _ ^ )(y _ ^ ' ( 2 ~ 22 ) 



i=l n=0 ii<i2<...<i n 

This formula allows us to obtain the following expression for the expectations 

1 1 



TV 



x - Mi y - M 2 J I VliV2 
at] r N "+ 1 -, 

- Z~ s / n*W%) A(X)A(r)eE-«E(-l)"E n F^lfe^J 

JV J i=l n=0 ii<...<i„+i fc=l V fe;vy y fc; 

1 N-l N n+1 

= w £ (- 1 ) n E e ^) E / II 4>«u)(*i)<t>T U )(yj)e XM II ^7— y ( 2 - 2 3) 

' n=0 ff.reSjv ii<i 2 <...<i„ + i j=l k=l V * fcAy ^ 

In eq. ()2-23|) we have used the normalized quasi-polynomial differentials defined in 1)1-71) . the fact that, with 
our normalizations for the Haar measure of U(N), the partition function is = (Nl) 2 Iljio 1 hj ( see ' ' j • 
and the identities 

N / N-l \ N 

a(x) n dp( Xk ) = n e <°) n ^mo**) ( 2 - 24 ) 

k=l \j=0 J crGSjv fc=l 

N I N-l \ JV 

A(y) n d^y fc ) = n \/m e e ( r ) n ^w(^) > ( 2 - 25 ) 

fc=l \ j=0 / re5jv fe=l 

which are obtained by replacing the monomials in the Vandermonde determinants by the biorthogonal polyno- 
mials of the same degree. 

2.2.1 Proof of formula <\TTE\) 

We are now in the position of proving formula ()1-16|) in a few strokes. Taking the coefficient of x r y s from the 
formal generating function in eq. (|2-2cS|) we obtain the expression 
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(Tr(M 1 r M 2 s )) Vi!V2 = 

min(r,s) N 



n=0 il<l2<— <«n+l o-.tS^jv i=l 

min(r,s) TV 



m E E E e (^)/ II **u)(*i)<iw)(vs)* iVi x ( 2 - 27 ) 

n=0 ii<l2<...<« n +l (t,t65jv j£{ii,...,i n +i} 



In this formula the notation xu-\ x means the sequence of variables x^, . . . ,Xi n+1 (similarly for the y's) and 
we have used the fact that the Schur polynomials S r as defined in (|2-18j) vanish if the number of variables is 
greater than r. Next, the orthogonality relations between the n 's and the VVi's in line Q2-27|) imply that the 
sum over the permutations a, r is restricted to those permutations such that 

o- = 7-077, 77 E S{h, ■ ■ ■ i n +i} (2-29) 

where S{i±, . . . i n+ \} denotes the group of permutation of the indices ik- The restriction on the indices 
il, . . . ,i n +i can be lifted because the following expression is permutation invariant in the label of those in- 
dices and when two such indices coincide the corresponding term vanishes due to the alternating form of the 
sum. This will produce an overcounting of a factor (n + 1)! which must be corrected: moreover we can relabel 
the variables of integration of the integral in line (|2-28j) from x,- lh , yi k to xt , Vk 

(Tr(M 1 r M 2 s )) Vi>V2 = 

min(r,s) , .„ ^ n+1 

= ^E E ^ E E^)/ eEr ^^ 

o-eSjv n =0 ' ii,t2,...,in+i i)£S„ + i k=l 

min(r,s) ,_-,x n „ n+1 

n=0 ^ U,»2,— i*n+l ^eSn+l k=l 

min(r,s) /_..%„ n+1 

= E T^+\)\ E E 6 ^)/ ^2 +1 ^5 r ^ [lj ... in+1] )5 s (y [1) ... jn+1] )n^(^)^ Cfe) (^)= (2-32) 

n=0 ^ ' ii,t2,...,»n+i i)6S„ + i fc=l 

min(r,s) ,_-,x n „ n+1 

E T^rW E E e ^)/ ^^^[i,...,^])^^^ 



n=0 U,«2,— !»n+l »7eSVi+i fc=l 

min(r,s) 



E T^Tv/ e E2+laik, ' k 5r(x [1> ... in+1] )5.(y [1 ,„. >n+1] )det 

n=0 ^"^^^ 



Ki2(xj,ye) 



(2-34) 

j,^<n+l 



where we have used the definition of the kernel Ki2(x,y) given in eq. (|1-19|) . This concludes the proof of eq. 
from which formula (jl-18|) follows immediately. 
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2.2.2 Proof of formula (Trill 

Resuming from eq. (|2-23|) and performing a relabelling of the x's and the y's allows to choose = k, and the 
sum over i\ < . . . < i n +\ becomes a combinatorial factor: 



Tr 



1 1 

x-M iy -M 2/ , u (j 



1 / AT \ f 1 

n=0 V / a,reS N j=l fe=l V ny y ' 

Now, the equations 1)1-51) and (|1-8|) imply that 
and 

/ ^ x ')^){yY V \ x _J ){y _ yl) = W„ U) , Tij) (2-37) 
where W is the N x N square matrix 4 : 

W := 7T — 1 . (2-38) 
nx -Qy - Pn 

Therefore t = or) where r/ is a permutation of the n + 1 first indices only, and we can write: 

n+l 

"7(i) 



Tr 



1 1 



x-M iy -M 2 J / Vi V2 



«=0 v y aeS N i;eS„ + i j=i 

E (n + l)l(JV -n-l)! E det W,,*! ■ ( 2 ~ 39 ) 



We note a,- = a(i) and we write: 



Tr 



1 1 \ \ °° (— 1 ^ n 

) / W .» = S (" + -"-')' j~L det ^^ »«• (2 " 40) 



x - Mi y - M 2 , 



v 1>V2 n=0 V*t w " ^W-^ 
where the sum over a n +2, • • • , ckjv disappears and brings a factor (iV — n — 1)!: 



Tr 



1 1 \\ ( — l) n 

x _Miy-M 2 )/ v v = E (n + l)! E det (^>%')i,;=l,.,n+l 

OO / -,\n n+l 



= e e ^n^w ( 2 - 41 ) 

n=0 V a 1 ^...^a n+ i ?7G5„+i j=l 

4 We remind the reader that this matrix should be properly understood as a formal power series in inverse powers of x and y, 
although an analytic definition could be given in terms of the biorthogonal polynomials. However this is unnecessary for the scope 
of the present paper. 
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and we can replace the sum over distinct a's by a sum over all a's: 



Tr 



1 



1 



x - Mi y - M 2 



Et^At E ^ E II *W 



n+1 



(n + 1) 



(2-42) 



Vi,V 2 n =0 K " ' ''r]£Sn+i ai,...,a„+i j=l 

If n is decomposed into a product of M{rj) cyclic permutations of lengths h(j]) + . . . + = n + 1, we have: 



Tr 



1 



1 



x - Mi y - M 2 



oo , s n Af(v) 

= Et^tv e m-v iM+i ™ iAv) 

V U V 2 n=0^ + i \eS n+1 j=l 

oo Af(n) 

-e^e ni-™- '"). 

n=0 v ; 'r;G5 n+ i j=l 

Now we use the following Lemma, which is a classical result in combinatorics: 



(2-43) 



Lemma 2.1 Let Q m be a function defined on the permutation group of m elements S m with the cluster property, i.e., 
such that if T) = r/i o ry 2 is a decomposition into disjoint permutations of m! and m" elements ( and hence m = m! + m") 
then 

G m (v) =Qm>{Vl)Gm»{V2) ■ (2-44) 

Under these circumstances we have the identity 



ex p I E E tj sm = 1 + E ^ E g "^) , 

I m=l ves m 



(2-45) 



\m=lat£C m / 

where C m denotes the set of all permutations of maximal length and has cardinality (m — 1)!. 



In other words this lemma says that if Q has the cluster property, then taking the logarithm of the RHS of 
eq. (|2-45p removes all "nonconnected" contributions and returns the "connected components" only. In view of 
Lemma 12. II let us define 

Gm(v) ■= II (2-46) 

3=1 



which has clearly the cluster property, and we have: 

1 1 



1 - ( Tr 



x -Mi y - Mi, 



)) =i + J2^iY, G ^ = ex p 

' I Vi,V2 m =l jj^s^ 



oo 

E E 



m=l nGC„ 



exp 



E 

m=l 



TrW r 



exp [Trln (1 N -W)] =det(l N -W) 



(2-47) 



This concludes our proof of formula (|1-11|) . 



3 Conclusions 

Formula l|l-ll|) is quite simple in spite of the long computations involved in its proof. It is tempting to imagine 
that also more complicated multi-correlators could be reduced to a computation involving the matrices P and 
Q, i.e., to biorthogonal polynomials. Quite clearly the possibility rests on having a manageable formula for 
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the generic multi-correlators of the Itzykson-Zuber-Harish-Chandra integral. There are indications that such a 
formula should be derivable: for example it is quite simple to obtain a formula for a correlator of entries with 
the same first index. Indeed starting from eq. (|2-7|) and proceeding in the same way we did in order to obtain 
eq. ()2-15|) one can easily prove the following formula 



det 



U(N) 



£,m=l...N 



a=l 



(3-1) 



A(X)A(Y)HH(x i -x ia )l[(y £ -y j ) 

a=l £jti a 



where 



{ym-yjT n^-^)-(^-%) n_i e n^-^)+- • -+(-i) n ^ 



a=l 



s=0 ^ a=l 



(3-2) 
(3-3) 



However the computation for non-equal first indices becomes quickly extremely complicated at least using the 
technique in [T7|. Nonetheless we hope that this first computation can break through the general belief that 
computations of correlators in multi-matrix models are not feasible in the finite-N regime due to the angular 
integrations. 

It should also be remarked that for polynomial potentials we could use the so-called "loop" equations to 
obtain information on other correlators. The result, however would be dependent on the specific form of the 
potentials and would not provide information on the HCIZ integral itself. 

Let us also mention that this calculation could be generalized to other random matrix models, in particular 
the complex matrix model, which presents many similarities with the 2-matrix model. The gaussian complex 
matrix model has attracted lot of attention in string theory. A particular case of the ADS/CFT, is the 
conjectured duality between string-theory in a pp-wave background, and BMN gauge theory. The gaussian 
complex-matrix model appears as an effective BMN theory in a particular limit, and the computation of mixed 
correlation functions is very important in that model [!][§]. In the gaussian complex matrix model, a formula is 
known for the 2-point mixed correlator (Tt x ^ m _} M ^ ), but little is known for other mixed correlation functions. 



A Equivalence of eq. ( l^HH and ( l^T5ft 



In this appendix we prove that the two formula? 
(|2-15|) and compute 



(\ u ji\ 2 )u{N) 



and (|2-15j) are equivalent. To this end we start from 



det 



(xi - Xi)(y m - yj) - 1 )e XiV 



'm=l...N 



A(X)A(Y) 



E <p)-- 



Y[( x t - x i)Y[(V£ - Vj) 



a{x)a{y)^ n ^u^-^n^-w) 



(A-l) 



fl ((x e - Xi)(y p{e) - Vj ) - l) = (A-2) 
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N 



t+i 1+3 



A w A ( y )p£ e(p) n^ : ^n^-^ 

1 - 2 (a;/ - Xi)(y p(e) - y 3 ) + E (x^ - Xi){y p ( tl ) - Vj){xt 2 - Xi)(y p{ e 2) -%■) + ... 
Let us consider the following subexpression from the above formula 



(A-3) 
(A-4) 

(A-5) 
(A-6) 

(A-7) 



_1\AT+1 



(-1) 

Y[(xe - Xi)Y[(y e - yj) 

t+i i+3 



N n 

i + ^2(-i) n E II( x 4-^)(yp(4)-%) 



n=l *i<* 2 < — <*n fc = l 



W 



(A-8) 



JV-l 

i + ^(-l) n E 

n=l e 1 <e 2 <-- 



n 

Y[(x tk - Xi)(y Ktk) - yj) 



k=i 



^<e 2 <...<i n ]J(xt - Xi) Y[(y p (e) - Vj) 

i+i e+j t k *i,p(t k w*k 



(A-9) 



(-1) 



N+l 



N-l 



v 



Y[(x<> -x^Yi (vt-yj) n =i «i<«2<-<^ n o^-^) n (vpw-vj) 



^{i,4,Vfc} pW0{i,p(4),Vfc} 



(A-10) 



(-1) 



JV-l 



Y[(xt - x^ Y[(y e - yj) 



N-n+l 



E 



n=l 



n - ^) n ~ 



<A-11) 



("I) 



iV+1 



N-l 



+ E(- 1 ) n+1 E 



i+i 1+3 ie{i k }je{p(i k )} k:ik ^ k: P (i k )+j 



KA-12) 



14 



N 



E(-D n+1 E 



i 



i iK...<i„ : Ylfa h -xi) JJ (y p ( ik )-yj) 



=(A-13) 



(**) 



E(-!) n+1 E- 



8=1 



n=l 



«l<...<«r7 



A(x il ,...,Xi n )A(y p(n) ,...,y p(in) ) 



=(A-14) 



det 



A? 



E(- 1 ) n+1 E 



1 \ / 1 



a;"" 1 



det 



71=1 



2l<...<«„ 



det 



/ 1 



1 \ 



y p{h) 

\°3P(ii) 



X 



det 



/ 1 

VM.ii) 



1 \ 

Vp(in+l) 



n-1 
Vp(i n +i) 

fijp(in + l) / 



-J 



1 \ 

Vp{in + l) 



(A-15) 



V^p(ii) "' Vp(i n+1 )J 



In the above chain of equality we have replaced (after the (*)) the upper limit of summation by N — 1 because 
in the case n = N there is certainly one £k equal to i (and one p(lk) equal to j) so that that term does not 
contribute. After (**) we have removed the condition on the multi-index because it is implicit in the sum of 
deltas that follows. Putting together eq. ()A-15|) with eq. (|A-7|) or (|2-15|) we obtain the desired proof of the 
equivalence of formula ()2-6j) with eq. (|2-7j) , thus also proving the former rigorously, which was not done in ^B] • 
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